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The chiral two-pion exchange nucleon-nucleon interaction has nowadays a rather firm conceptual 
lO ■ basis, but depends on low-energy constants, which may be extracted from fits to data. In order to 

' facilitate this kind of application, we present here a parametrized version of our relativistic expansion 

of this component of the force to 0{q'^), performed recently. 
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I. INTRODUCTION 



Nuclear interactions are strongly dominated by the quarks u and d and can be accomodated into 
a two-flavor version of QCD. The masses of these quarks are small in the GeV scale and one is 
^ ' not far from the massless lagrangian limit, in which QCD is invariant under chiral SU{2) x 5*17(2) 

^ : transformations. For this reason, low energy hadronic processes can be reliably described by means 

of effective lagrangians that are symmetric under the Poincare and isospin groups and incorporate 
approximate chiral symmetry, realized in the Nambu-Goldstone mode. 

The one-pion exchange A'^A'^ potential {OPEP) is simple, has been well established long ago, and 
dominates completely partial waves with orbital momentum L > 5. The two-pion exchange potential 
, (TPEP), on the other hand, is rahter complex and has become free of important umbiguities only 

in the 1990s, after the systematic use of chiral symmetry in its theoretical description [1-11]. 

Chiral perturbation theory is based on the existence of a characteristic scale q, set by both pion 
four-momenta and nucleon three-momenta, such that q < 1 GeV. Due to this technique, nowadays 
^ ' one understands rather well the internal hierarchies of the A'^A^ potential in terms of chiral layers. 

'~i \ Leading terms of the chiral TPEP are of order Oijf) and expansions which go up to 0{q'^) are 

P5 . already available. One of them was produced recently by our group [10,11]. We departed from a 

relativistic lagrangian and evaluated the relevant Feynman diagrams covariantly, without resorting 
to heavy baryon approximations. The so obtained T-matrix was then transformed into a potential, 
■ expressed in terms of covariant loop integrals and observable parameters. Without loss of generality, 

5h ' one may choose these parameters to be either the subthreshold coefficients extracted from nN 

scattering [12] or the low-energy contants (LECs) present in the effective lagrangian. These two 
possibilities are formally equivalent, but suitded to slightly different physical purposes. The former 
choice yields a closed prediction for the potential, whereas the latter gives rise to an open theoretical 
structure which may be used to fit A'^A' data. 

Nowadays, a rather pressing issue is to determine the extent that nature backs this picture. The 
comparison of chiral predictions with empirical phase shifts is not straightforward, since existing 
theoretical potentials are not reliable for distances smaller than 1 fm [11]. Three complementary 
possibilities are available for overcoming this problem. The most direct one is to use peripheral 
waves and rely on those windows in angular momenta and energies for which the centrifugal barrier 
is effective in suppressing the interaction at short distances. When this happens, the Born approx- 
imation can be used and one does not need to know the potential close to the origin. However, 
this kind of test is not very stringent, since peripheral waves are small and uncertainties are large 
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[4,5,13]. Inclusion of more important waves requires the use of dynamical equations with regularized 
potentials [8]. This regularization brings necessarily extra parameters into the problem which are 
not constrained by chiral symmetry. One is then faced with the problem of disentangling from ones 
results those windows which do indeed test the symmetry. The third alternative, already employed 
by the Nijmegen group [7], is to assume that the theoretical potential determines correctly the in- 
teraction in a spatial window ranging from a radius R onwards and then use it as an input in phase 
shift analiscs. This procedure has already proved to be useful in the case of the OPEP, in the 
clastic regime and for waves with L > 5. Its extension to the TPEP becomes not trivial due to 
the fact that the range of reliability of the theoretical potential depends on the chiral order one is 
working at. 

The research on the TPEP performed in the last decade has set its conceptual foundations on 
a rather solid basis, comparable to that of the OPEP in the late sixties. On the other hand, the 
TPEP depends on several LECs, which must be extracted from eihter ttN scattering data or direct 
fits of NN phase shifts. In general, this last kind of procedure tends to be computationally heavy, 
for theoretical results are usually given as cumbersome expressions. In order to make applications 
easier, in this work we present a parametrized version of our 0(g*) relativistic configuration space 
TPEP, which is numerically accurate for distances larger than 1 fm. It is based on the theoretical 
expressions derived in Ref. [11] and reproduced in appendix A, as functions of the LECs. 



II. PARAMETRIZED POTENTIAL 



The configuration space potential has the isospin structure 



V{r) = V+{r) + T'^^^-T^^^ V-{r) 



with 



V^ir) = V^+ V^s ^Ls + V^aT + Qss + fig , 



and Qls = L-{ct'^'> +(t''^>)/2, Qt = 3ct^' -r cr'^' -r - ct^'> -a^^' , Qss = (t^' -cr^' ■ The iorm of 
the operator fig in configuration space is highly non-local and can be found in Ref. [14]. 

In Refs. [10,11] we have presented a 0{q^) relativistic expansion of the TPEP, which is repro- 
duced, in an alternative form, in appendix A. The configuration space potential is written in terms 
of numerical coefficients which multiply dimensionless functions arising form the Fourier transforms 
of Feynman loop integrals. The former are combinations of external parameters representing the 
pion and rmcleon masses, /^ and m, respectively, the pion decay constant /,r, the axial coupling 
constant qa, and the LECs Ci and di. The latter are denoted by Zi and depend on just n, m, /„, 
and X = fj,r. We keep the external quantities as free and parametrize the function Zi = {Fi, Gi) as 



Zi = 



(47r)-'/2 



7i X 



The coefficients 7" corresponding to the various cases are given in the tables at the end of this 
section. This parametrization is more than 1% accurate in the range 0.8 fm < r < 10 fm. 
Using the definition a = /x/m, the profile functions are written as 
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The parametrized profile functions given above depend explicitly on four well known quantities, 
namely m, /i, qa, /tt, and on the less known LECs Ci and di. Therefore, the latter may be extracted 
from fits to data. When doing this, however, one has to bear in mind that, as discused in Ref. [11], 
the influence of the LECs over the profile functions is rather uneven. Indeed, their influence over 
Vp", V^g, VjT, and Vgg is rahter strong, but barely perceptible in VJ" , V^, and V^^. 
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APPENDIX A: THEORETICAL POTENTIAL 

The 0(g*) relativistic expansion of the TPEP produced in Refs. [10,11] was based on the eval- 
uation of three familes of diagrams^. The first of them involves only pion and nucleon degrees of 
freedom into single loops and corresponds to the minimal realization of chiral symmetry [3]. It 
includes the subtraction of the iterated OPEP and yields the terms in the profile functions given 
below which are proportional to just QaI ft, Qa/ ft or 1//*. Terms proportional to 1//^, on the 
other hand, come from two-loop processes, either in the form of f-channel contributions from the 
second family or s and w-channel terms embodied in the subthreshold coefficients of the third fam- 
ily. Finally, the third group of diagrams includes chiral corrections associated with other degrees of 
freedom, hidden within the LECs Cf and di, and gives rise to contributions which are proportional 
to either {LEG)/ ft or (LEC)'^ /ft- In Ref. [11] we have expressed this last class of results in terms 
of the nN subthreshold coefficients. As these can be easily translated into LECs, in the present 
work we write the potential in terms of these constants, which appear directly into the effective 
lagrangians. The following expressions correspond to the updated version as described in Ref. [13]. 

The radial components of the potential are expressed in terms of the following profile functions'^ 
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where r"*" = 3 and t~ = 2. 

Defining t as the Laplacian operator acting on the variable x = ij,r we write our profile functions 

as 
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^Please see section II of Ref. [11], for a detailed discussion of the meaning and dynamical contents of these families of diagrams, 
which are given in its Fig. 2. 
^Please see Eqs. (3.4)-(3.8) of Ref. [11]. Note that in Eq. (3.8) a multiplication factor of is missing. 
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The dimensionless functions Si {x) carry the spatial dependence of the potential and are given by 
Eqs. (3.16)-(3.23) of Ref. [11]. 
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